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Abstract 

We propose an adaptive transmission technique for free space optical (FSO) systems, operating in 
atmospheric turbulence and employing subcarrier phase shift keying (S-PSK) intensity modulation. Ex- 
ploiting the constant envelope characteristics of S-PSK, the proposed technique offers efficient utilization 
of the FSO channel capacity by adapting the modulation order of S-PSK, according to the instantaneous 
state of turbulence induced fading and a pre-defined bit error rate (BER) requirement. Novel expressions 
for the spectral efficiency and average BER of the proposed adaptive FSO system are presented and 
performance investigations under various turbulence conditions and target BER requirements are carried 
out. Numerical results indicate that significant spectral efficiency gains are offered without increasing 
the transmitted average optical power or sacrificing BER requirements, in moderate-to-strong turbulence 
conditions. Furthermore, the proposed variable rate transmission technique is applied to multiple in- 
put multiple output (MIMO) FSO systems, providing additional improvement in the achieved spectral 
efficiency as the number of the transmit and/or receive apertures increases. 
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I. Introduction 

Free Space Optical (FSO) communication is a wireless technology, which has recently attracted 
considerable interest within the research community, since it can be advantageous for a variety of 
applications HI- Q- However, despite its significant advantages, the widespread deployment of FSO 
systems is limited by their high vulnerability to adverse atmospheric conditions |4). Even in a clear sky, 
due to inhomogeneities in temperature and pressure changes, the refractive index of the atmosphere varies 
stochastically and results in atmospheric turbulence. This causes rapid fluctuations at the intensity of the 
received optical signal, known as turbulence-induced fading, that severely affect the reliability and/or 
communication rate provided by the FSO link. 

Over the last years, several fading mitigation techniques have been proposed for deployment in 
FSO links to combat the degrading effects of atmospheric turbulence. Error control coding (ECC) in 
conjunction with interleaving has been investigated in [5] and [6]. Although this technique is known 
in the Radio Frequency (RF) literature to provide an effective time-diversity solution to rapidly-varying 
fading channels, its practical use in FSO links is limited due to the large-size interleave^ required to 
achieve the promising coding gains theoretically available. Maximum Likelihood Sequence Detection 
(MLSD), which has been proposed in [7], efficiently exploits the channel's temporal characteristics; 
however it suffers from extreme computational complexity and therefore in practice, only suboptimal 
MLSD solutions can be employed HI- ifTOll . Of particular interest is the application of spatial diversity to 
FSO systems, i.e., transmission and/or reception through multiple apertures, since significant performance 
gains are offered by taking advantage the additional degrees of freedom in the spatial dimension [11]- 
lfT2ll . Nevertheless, increasing the number of apertures, increases the cost and the overall physical size 
of the FSO systems. 

Another promising solution is the employment of adaptive transmission, a well known technique 
employed in RF systems 11131 . By varying basic transmission parameters according to the channel's 
fading intensity, adaptive transmission takes advantage of the time varying nature of turbulence, allowing 
higher data rates to be transmitted under favorable turbulence conditions. Thus, the spectral efficiency 
of the FSO link can be improved without wasting additional optical power or sacrificing performance 
requirements. The concept of adaptive transmission was first introduced in the context of FSO systems in 
lfl4ll . where an adaptive scheme that varied the period of the transmitted binary Pulse Position modulated 

'For the signalling rates of interest (typically of the the order of Gbps), FSO channels exhibit slow fading, since the correlation 
time of turbulence is of the order of 10~ 3 to 10~ 2 seconds. 
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(BPPM) symbol was studied. Since then, various adaptive FSO systems have been proposed. In lfT5l . a 
variable rate FSO system employing adaptive Turbo-based coding schemes in conjunction with On-Off 
keying modulation was investigated, while in |fT6ll an adaptive power scheme was suggested for reducing 
the average power consumption in constant rate optical satellite-to-earth links. Recently, in iTTTI . an 
adaptive transmission scheme that varied both the power and the modulation order of a FSO system with 
pulse amplitude modulation (PAM), has been studied. 

In this work, we propose an adaptive modulation scheme for FSO systems operating in turbulence, 
using an alternative type of modulation; subcarrier phase shift keying intensity modulation (S-PSK). 
S-PSK refers to the transmission of PSK modulated RF signals, after being properly biasedj, through 
intensity modulation direct detection (IM/DD) optical systems and its employment can be advantageous, 
since: 

• in the presence of turbulence, it offers increased demodulation performance compared to PAM 
signalling HI- 1H, 

• due to its constant envelope characteristic, both the average and peak optical power are constant in 
every symbol transmitted, and, 

• it allows RF signals to be directly transmitted through FSO links providing protocol transparency 
in heterogeneous wireless networks 11211 - ll23ll . 

Taking into consideration that the bias signal required by S-PSK in order to satisfy the non-negativity 
requirement is independent of the modulation order, we present a novel variable rate transmission strategy 
that is implemented through the modification of the modulation order of S-PSK, according to the 
instantaneous turbulence induced fading and a pre-defined value of Bit Error Rate (BER). The performance 
of the proposed variable rate transmission scheme is investigated, in terms of spectral efficiency and 
BER, under different degrees of turbulence strength, and is further compared to non-adaptive modulation 
and the upper capacity bound provided by |[24l . Moreover, an application to multiple input multiple 
output (MIMO) FSO systems employing equal gain combining (EGC) at the receiver is provided and the 
performance of the presented transmission policy is evaluated for various MIMO deployments. 

The remainder of the paper is organized as follows. In Section [TTJ the non-adaptive S-PSK FSO system 
model is described and its performance in the presence of turbulence induced fading is investigated. In 
section [TTlJ the adaptive S-PSK strategy is presented, deriving expressions for its spectral efficiency and 

2 Since optical intensity must satisfy the non-negativity constraint, a proper DC bias must be added to the RF electrical signal 
in order to prevent clipping and distortion in the optical domain. 
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BER performance, and an application to MIMO FSO systems is further provided. Section HVl discusses 
some numerical results and useful concluding remarks are drawn in section [V] 

Notations: E{-} denotes statistical expectation; N (/U,cr 2 ) denotes Gaussian distribution with mean \i 
and variance a 2 . 

II. Non-Adaptive Subcarrier PSK Intensity Modulation 

We consider an IM/DD FSO system which uses a subcarrier signal for the modulation of the optical 
carrier's intensity and operates over the atmospheric turbulence induced fading channel. 

A. System and Channel Model 

On the transmitter end, we assume that the RF subcarrier signal is modulated by the data sequence 
using PSK. Moreover a proper DC bias is added in order to ensure that the transmitted waveform always 
satisfies the non-negativity input constraint. Hence, the transmitted optical power can be expressed as 

P t (t) =P[l + fiS (t)] (1) 

where P is the average transmitted optical power and \i is the modulation index (0 < \i < 1) which 
ensures that the laser operates in its linear region and avoids over-modulation induced clipping. Further, 
s (t) is the output of the electrical PSK modulator which can be written as 

s (t) = J2 9 (t-kT) cos (2nf c t + cj) k ) (2) 
k 

where f c is the frequency of the RF subcarrier signal, T is the symbol's period, g (i) is the shaping pulse, 
4>k £ [0, (M — 1) jjl is the phase of the fcth transmitted symbol and M is the modulation order. 

On the receiver's end, the optical power which is incident on the photodetector is converted into an 
electrical signal through direct detection. We assume operation in the high signal-to-noise ratio (SNR) 
regime where the shot noise caused by ambient light is dominant and therefore Gaussian noise model is 
used as a good approximation of the Poisson photon counting detection model Q. 

After removing the DC bias and demodulating through a standard RF PSK demodulator, the electrical 
signal sampled during the kth symbol interval, which is obtained at the output of the receiver, is expressed 

as 

r [k] = nn^J^-PI [k] s [k] + n [k] (3) 

where i] corresponds to the receiver's optical-to-electrical efficiency, s [k] = cos 4>k — j sinc^, E g is the 
energy of the shaping pulse and n [k] is the zero mean circularly symmetric complex Gaussian noise 
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component with E{n[/c]n* [k]} = 2a^ = N a . Furthermore, I[k] represents turbulence-induced fading 
coefficient during the kth symbol interval and is given by 

/ [k] = I exp (2x [k]) (4) 

where I a denotes the signal light intensity without turbulence and x [k] is a normally distributed random 
variable with mean m x and variance a x , i.e. j x n.\ (x) = N (m x ,a x ). Hence I [k] follows a lognormal 
distribution with probability density function (PDF) provided by 



( < <^ ^ 



\ 



In Ij-J ~2m x 
8cj2 



(5) 



To ensure that the fading does not attenuate or amplify the average power, we normalize the fading 
coefficients such that E j | = 1. Doing so requires the choice of m x = —a x ifTTTl . Moreover, 
without loss of generality, it is assumed that I = 1. 

Atmospheric turbulence results in a very slowly-varying fading in FSO systems. For the signalling rates 
of interest ranging from hundreds to thousands of Mbps |[25l . the fading coefficient can be considered 
as constant over hundred of thousand or millions of consecutive symbols, since the coherence time of 
the channel is about l-100ms |[26l . Hence, it is assumed that turbulence induced fading remains constant 
over a block of K symbols (block fading channel), and therefore we drop the time index k, i.e. 

I = I[k], k = l,...K (6) 

It should be noted that in the analysis that follows, it is further assumed that the information message is 
long enough to reveal the long-term ergodic properties of the turbulence process. 
The instantaneous electrical SNR is defined as 

AW 

7= N n (7) 



while the average electrical SNR is given by 



with E s = ^. 



B. BER Performance 

The BER performance of the FSO system under consideration depends on the statistics of atmospheric 
turbulence and the modulation order. 
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When M = 2 (BPSK), the conditioned on the fading coefficient, I, BER is given by [20] 

P b (2,I) = q(i^) (9) 



while for M > 2 the following approximation can be used IT271 

Pb ^ I) -^i Q ( I ^^ : ¥i) (10) 

where Q (•) is the Gaussian Q-function defined as Q (x) = |™e 3"(ix. Hence, the average BER 
will be obtained by averaging over the turbulence PDF, i.e. 

POO 

P b (M) = / P b (M,I)f T (I)dI (11) 
Jo 

which can be efficiently evaluated using the Gauss-Hermitte quadrature formula ifTTTl . 

C. Capacity Upper Bound 

Using the trigonometric moment space method, an upper bound for the capacity of optical intensity 
channels when multiple subcarrier modulation is employed, has been derived in 11241 . By applying these 
results to the FSO system under consideration (one subcarrier), the conditioned on the fading coefficient, 
/, capacity can be upper bounded by 

T / u 2 ^ 2 P 2 E,T 2 \ 1 

(12) 



, , l x 2 i 1 2 P 2 E g I 2 , 
log 2 7T + log 2 ' 9 +o{a n ) 



W 

Cu P (I) = y ^ e Ar o 

where W denotes the electrical bandwidth and o (a n ) represents the capacity residue which vanishes 
exponentially as a n — > 0. Hence at high values of electrical SNR, (fT2l) can be approximated by 

C up (I)*W\og 2 (?pJ. (13) 

The unconditional approximative capacity upper bound will be obtained by averaging (fT"3l) over the 
fading distribution, i.e. 

c up ~yJ °° log2 (^f ) fl (/) dI > (14) 

which, according to the Appendix, can be written in closed form as 

E ( log2 (1) _ M 

2 V \e) ln2 
and will be used as a benchmark in the analysis that follows. 

III. Adaptive Modulation Strategy 

In this section we introduce an adaptive transmission strategy that improves the spectral efficiency of S- 
PSK FSO systems, without increasing the transmitted average optical power or sacrificing the performance 
requirements. 



C up n — (log 2 (^) -^) (15) 
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A. Mode of Operation 

By inserting pilot symbols at the beginning of a block of symbols^, the receiver accurately estimates 
the instantaneous channel's fading state, /, which is experienced by the remaining symbols of the block. 
Based on this estimation, a decision device at the receiver selects the modulation order to be used for 
transmitting the non-pilot symbols of the block, configures the electrical demodulator accordingly and 
informs the adaptive PSK transmitter about that decision via a reliable RF feed back path. 

The objective of the above described transmission technique is to maximize the number of bits 
transmitted per symbol interval, by using the largest possible modulation order under the target BER 
requirement P a . Hence the problem is formulated as 

max logo M 

(16) 

s.t. P b (M,I)<P 

In practice, the modulation order will be selected from N available ones, i.e. {M±, M2, Mn}, de- 
pending on the values of I and P a . Specifically, the range of the values of the fading term is divided in 
(N + 1) regions and each region is associated with the modulation order, Mj, according to the rule 

M = Mj = 2 j if Ij <I < I j+ i, j = 1, ...N (17) 

The region boundaries {Ij} are set to the required values of turbulence-induced fading required to achieve 
the target P a . Hence, according to (O and ( flOl ), they are obtained by 

h = \[^=Q~ l ( p o) » (18) 

I j = ^flQ- 1 (^Po),j = 2,...N (19) 



' ^ V 27 



and 



In+i = R, (20) 

where R — > 00 and (•) denotes the inverse Q-f unction, which is a standard built-in function in most 
of the well-known mathematical software packages. It should be noted that in the case of / < I\, the 
proposed transmission technique stops transmission. 

3 Taking into consideration the length of a transmitting block of symbols, the insertion of pilot symbols will not cause significant 
overhead. 
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B. Performance Evaluation 

1 ) Achievable Spectral Efficiency: The achievable spectral efficiency is defined as the information rate 
transmitted in a given bandwidth, and for the communication system under consideration is given b)Q 

5=- = ^ (21) 

with C representing the capacity measured in bit/s and n the average number of transmitted bits. 
The average number of transmitted bits in the adaptive S-PSK scheme is obtained by 

N+l 

™ = ^2 a 3 l0g 2 M 3 ( 2T > 

where 

aj = Pr {Ij <I< = [ 3+1 fj (I) dl. (23) 

Taking into consideration that the CDF of the LN fading distribution is given by 

eq. d23l can be equivalently written as 

Oj = Fj (J i+ i) - Fj (Ij) = Q ( Xj ) - Q (x j+ i) (25) 

with xj = ln ^2^ 2Ul . Hence, the achievable spectral efficiency can be written as 

Ef=i ajlog 2 Mj E?=iQ(*j)-(N + l)Q(x N+1 ) 

S= — 2 = 2 ' (26) 

which is simplified to 

S = 3 ^ (27) 

since Q (a;jv+i) — > 0, according to (l20l ). 

2) Average Bit Error Rate: The average BER of the proposed variable rate FSO system can be 
calculated as the ratio of the average number of bits in error over the total number of transmitted bits 

|[T3l . The average number of bits in error can be obtained by 

N+l 

Uerr = ( Pb h l0 §2 M 3 ( 28 ) 

0=1 

4 Note that since the subcarrier PSK modulation requires twice the bandwidth than PAM signalling, the average number of 
bits transmitted in a symbol's interval will be divided by two. 
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where 

(Pb), = £ Pb (Mj , I) fi (I) dl. (29) 
Hence the average BER is given by 

Pb = (30) 

n 

C. Application to MIMO FSO systems 

Consider a Multiple Input Multiple Output (MIMO) FSO system where the information signal is 
transmitted via F apertures and received by L apertures. For the MIMO system under consideration, it 
is assumed that the information bits are modulated using S-PSK and transmitted through the F apertures 
using repetition coding |[28l . Thus, the received block of symbols at the Zth receive aperture is given by 

n [k] = WP ^ S [k] J2 JC/O + in[k],k= 1, ..K (31) 
/=i 

where I^ l > denotes the fading coefficient that models the atmospheric turbulence through the opti- 
cal channel between the /th transmit and the /th receive aperture and n [k] represents AWGN with 
E {n [k] n* [k]} = 2o\. After equal gain combining (EGC) the optical signals from the L receive apertures, 
the output of the receiver will be obtained as 

r [k] = £ n [k] = ^>/£°w 22 j2 lUl) + n w • w 
i=i f=i i=i 

Note that the factor F that appears in (|3"TT ) and (l32l . is included in order to ensure that the total transmit 
power is the same with that of a system with no transmit diversity, while the factor L ensures that the 
sum of the L receive aperture areas is the same with the aperture area of a system with no receive 
diversity. Moreover, the statistics of the fading coefficients of the underlying FSO links are considered to 
be statistically independent; an assumption which is realistic by placing the transmitter and the receiver 
apertures just a few centimeters apart [26]. 

The variable rate subcarrier PSK transmission scheme can be directly applied to the MIMO configu- 
ration with the decision on the modulation order to be based on 

FL • <33) 

Hence, after determining the region boundaries for the target P requirement, using Eqs. (fT8l)-(|20l>, the 
optimum modulation order will be selected from the N available ones depending on the value of It, i.e. 

M = Mj if Ij <I T < I j+ i, j = 1, ...AT (34) 
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Similarly to the Single Input Single Output (SISO) configuration, the achievable spectral efficiency of 

the adaptive MIMO FSO system will be obtained by 

Ef-V bj log 2 Mj 
S = ^°~ x 3 3 (35) 

2 

where 

bj = Pr {Ij <I T < = [ 3+1 f lT (I) dl (36) 

J I, 

and fj T (I) is the the PDF of It, which can be efficiently approximated by the lognormal distribution 
EE EH, i.e. 

I T = exp (0 (37) 

with (£) = N (rn£,crg\, cr| = In (l + ^—prT^ anc ^ m 5 = ~\ a \- Hence, Eq. d36l ) can be equivalently 
written as 

bj = Q {Vj) ~ Q (Vj+i) (38) 
with yj = — 2 and d35l) is reduced to 

5 = SiM. (39) 
Furthermore, the average BER of the adaptive MIMO FSO system will be obtained by 



h = i^+i . - I : (40) 



EgT log 2 M j 
Ef = V 6i log 2 M. 
where 

(P 6 ), = / P b (M„ I) / /t (/) dl. (41) 



IV. Results & Discussion 

In this section, we present numerical results for the performance of the adaptive S-PSK scheme in 
various turbulence conditions and for different target BERs. We further apply this transmission policy to 
different MIMO deployments. 

Figs. [T][3] depict the spectral efficiency of the adaptive subcarrier PSK transmission scheme at different 
degrees of turbulence strength, i.e. a x = 0.1, a x = 0.3 and a x = 0.5, and when a SISO FSO system is 
considered. Specifically, numerical results obtained by d27T ) for two different target BER requirements, 
P = 10~ 2 and P = 10~ 3 , are plotted as a function of the average electrical SNR, along with the upper 
bound provided by ([TBI ) and the spectral efficiency of the non-adaptive subcarrier BPSK (M = 2). The 
latter is found by determining the value of the average electrical SNR for which the BER performance 
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of the non-adaptive BPSK, as given by (TTTT >. equals P a . It is obvious from the figures that the spectral 
efficiency of the adaptive transmission scheme increases and comes closer the upper bound by increasing 
the target P . Moreover, when compared to the non-adaptive BPSK, it is observed that adaptive PSK 
offers large spectral efficiency gains (14dB when P Q = 10 -3 ) at strong turbulence conditions (a x = 0.5); 
however, these gains are reduced as a x reduces. For low turbulence (a x = 0.1), it is observed that 
non-adaptive BPSK reached its maximum spectral efficiency (S = 0.5) at lower values of SNR than the 
proposed adaptive scheme, indicating that in these conditions it is more effective to modify the modulation 
order based on the value of the average SNR rather than the instant value of fading intensity. Hence, 
the proposed adaptive PSK technique, which is based on the estimation of the instantaneous value of 
the channel's fading intensity, can be considered particularly effective only in the moderate to strong 
turbulence regime. 

Figs. |4jl6] illustrate the average BER performance of the adaptive transmission technique for the same 
target BER requirements and turbulence conditions. It is clearly depicted that the average BER of the 
adaptive system is lower than the target P Q in all cases examined, satisfying the basic design requirement 
of (fT6l ). Moreover, it can be easily observed that the performance of the adaptive system approaches 
the performance of the non-adaptive system with the largest modulation order, at high values of average 
SNR; this was expected, since in this SNR regime the adaptive scheme chooses to transmit with the 
largest available modulation order. 

Finally, Fig. [7] depicts numerical results for the spectral efficiency of various MIMO FSO systems, 
when the adaptive transmission technique with target P Q = 10~ 3 and N = 3 available modulation orders 
is applied and moderate turbulence conditions are considered. As it is clearly illustrated in the figure, 
the increase of the number of transmit and/or receive apertures improves the performance of the adaptive 
transmission scheme, increasing the achievable spectral efficiency. However, this does not happen at low 
values of average SNR (less than 8dB), which may seem surprising at first but can be explained by the 
following argument. At the low average SNR regime, the region boundaries {Ij} correspond to values 
higher than the unity. Hence, as the number of transmit and/or receive apertures increases, the variance 
o"! decreases, resulting in higher values for the parameters {yj} and decreased spectral efficiency. As 
the average SNR increases, most of the region boundaries {Ij} take values lower than unity and, as a 
concequence, the increase in the number of apertures results in lower values for the parameters {yj} and 
higher spectral efficiency. 
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V. Conclusions 

We have presented a novel adaptive transmission technique for FSO systems operating in atmospheric 
turbulence and employing S-PSK intensity modulation. The described technique was implemented through 
the modification of the modulation order of S-PSK according to the instantaneous fading state and a 
pre-defined BER requirement. Novel expressions for the spectral efficiency and average BER of the 
adaptive FSO system were derived and investigations over various turbulence conditions and target BER 
requirements were performed. Numerical results indicated that adaptive transmission offers significant 
spectral efficiency gains, compared to the non-adaptive modulation, at the moderate-to-strong turbulence 
regime (14dB at S = 0.5, when P Q = 1CT 3 and a x = 0.5); however, it was observed that in low 
turbulence, it is more efficient to perform adaptation based on the average electrical SNR, instead of the 
instantaneous fading state. Furthermore, the proposed technique was applied at MIMO FSO systems and 
additional improvement in the achieved spectral efficiency was observed at the high SNR regime, as the 
number of the transmit and/or receive apertures increased. 

Appendix 

This appendix provides a closed-form expression for (fT4l) . Using the pdf of turbulence induced fading, 
(fT4l) can be written as 

C up ^—I log 2 (l-)— ^==exp(-A — T J^)dI. (42) 



2 Jo e J 2/^2^1 "V 8a * 

To simplify (l42l . we substitute In / by y and hence 



where 



and 



C up ^K x +K 2 (43) 



^ilo,P)re X p(-^* (44) 



J-oc V 8a - 



W f 00 ( {y + 2al) 2 \ 



X 



Using HI Eq. (3.321/3)], (@4]> is reduced to 



2 oz Ve 

while, using (29l Eq. (3.461/2)], (g5j is reduced to 

2Wa 



K2 = - z T ^f. (47) 
In 2 



Hence, by substituting in (l43l) . Eq. (fT31) is obtained. 
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Fig. 1. Spectral efficiency of the adaptive subcarrier PSK scheme when a x = 0.1. 
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Fig. 2. Spectral efficiency of the adaptive subcarrier PSK scheme when a x = 0.3. 
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Fig. 3. Spectral efficiency of the adaptive subcarrier PSK scheme when a x = 0.5. 
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Fig. 4. Average BER of the adaptive subcarrier PSK scheme when a x = 0.1. 
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Fig. 5. Average BER of the adaptive subcarrier PSK scheme when a x = 0.3. 
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Fig. 6. Average BER of the adaptive subcarrier PSK scheme when a x = 0.5. 
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